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Tidal disruption has a dramatic impact on the outcome of neutron star-black hole mergers. The 
phenomenology of these systems can be divided in three classes: nondisrnptive, mildly disrnptive or 
disruptive. The cutoff frequency of the gravitational radiation produced during the merger (which 
is potentially measurable by interferometric detectors) is very different in each regime, and when 
the merger is disuptive it carries information on the neutron star equation of state. Here we use 
semianalytical tools to derive a formula for the critical binary mass ratio Q = Mbh/Mns below which 
mergers are disruptive as a function of the stellar compactness C = M-ms/Rns and the dimensionless 
black hole spin y. We then employ a new gravitational waveform amplitude model, calibrated to 
134 general relativistic numerical simulations of binaries with black hole spin (anti-)aligned with the 
orbital angular momentum, to obtain a fit to the gravitational-wave cutoff frequency in the disruptive 
regime as a function of C, Q and y. Our findings are important to build gravitational wave template 
banks, to determine whether neutron star-black hole mergers can emit electromagnetic radiation 
(thus helping multimessenger searches), and to improve event rate calculations for these systems. 

PACS numbers: 04.25.dk, 97.60.Jd, 97.60.Lf, 04.30.-w 


I. Introduction. The merger of black holes (BHs) and 
neutron stars (NSs) is one of the most violent events 
in the Universe. Coalescing NS-BH systems are among 
the leading candidate sources for upcoming interferomet¬ 
ric gravitational-wave (GW) detectors such as the Ad¬ 
vanced Laser Interferometer Gravitational-Wave Obser¬ 
vatory (LIGO) [1, 2], Virgo [3], the Kamioka Gravita¬ 
tional wave detector (KAGRA) [4, 5], and LIGO-India 
[6]. GW observations of NS-BH mergers may provide 
information on the NS equation of state (EOS) [7] and 
on the underlying theory of gravity [8]. NS-BH binaries 
are also short-hard gamma-ray burst (SGRB) progenitor 
candidates [9]. If the NS is tidally disrupted during the 
merger, a hot disk with mass > O.OIMq may form around 
the spinning remnant BH. A scenario where the BH-disk 
system launches a relativistic jet by releasing its grav¬ 
itational energy via neutrino or electromagnetic (EM) 
radiation on a time scale 2 s explains the duration, 
energetics, and estimated event rates of SGRBs [10, 11]. 
During their merger, NS-BH binaries may also emit EM 
radiation isotropically — as opposed to beamed SGRB 
emission — when they eject unbound material. This can 
be as massive as ~ O.IMq and have subrelativistic ve¬ 
locities of ^ 0.2-0.3c [12], producing EM counterparts in 
the form of macronov®/kilonov*, powered by decay heat 
of unstable r-process elements and by nonthermal radi¬ 
ation from electrons accelerated at blast waves between 
the merger ejecta and the interstellar medium [13-18]. 


The features of GW emission from NS-BH binaries, 
as well as the plausibility of these systems being cen¬ 
tral SGRB engines and sites for EM radiation emission 
in general, depend crucially on whether or not the NS is 
tidally disrupted. Only numerical simulations can assess 
this. Fortunately, with the enormous progress made over 
the last decade, numerical relativity has provided a clear 
picture of NS-BH GW emission and shedded light on the 
processes leading to disk formation and mass ejection. 
Most of the GW emission occurs prior to the NS tidal 
disruption, if this happens at all, and before significant 
thermal effects take place. Further, magnetic fields ap¬ 
pear to barely affect the GW signal. These circumstances 
imply that an ideal fluid-dynamics treatment with a cold 
EOS is appropriate to simulate the dynamical regime of 
interest for GW signals. A notable feature of these sig¬ 
nals is that the cutoff frequency at which their amplitude 
damps due to the NS tidal disruption depends on the NS 
EOS. Hence, the cutoff frequency encodes information 
on the EOS itself, particularly when this is stiff [19-22]. 

Numerical simulations of compact binary mergers are 
still very resource intensive, which makes semianalyti¬ 
cal models very valuable. Simulations and models are 
most advanced for BH-BH systems. Waveform models 
belong to two main classes: Fourier-domain phenomeno¬ 
logical inspiral-merger-ringdown (IMR) models based on 
a post-Newtonian (PN) description of the early inspi¬ 
ral stage (“PhenomA” [24-26], “PhenomB” [27], “Phe- 
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FIG. 1. The cutoff frequency /cut, as defined in Eq. (1), computed with our NS-BH GW amplitude model [23]. Each panel label 
specifies the NS EOS and the BH spin parameter y used. The contour lines report /cut in Hz and have a spacing of 250 Hz. 
The thick, green, continuous line is the location where the mass of the torus remnant Mb,torus vanishes. The two dashed lines 
in each panel divide the plane in three regions: a top-right region in which NS-BH coalescences are nondisruptive (/tide > /rd 
and Mb,torus = 0), a bottom-left one in which they are disruptive (/tide < /rd and Mb,torus > 0), and a middle region in which 
mildly disruptive coalescences occur (/tide < /rd and Mb,torus = 0, or /tide > /rd and Mb,torus > 0). 


nomC” [28], and “PhenomP” [29]) and effective-one- 
body models tuned to BH-BH simulations [30-36]. NS- 
BH waveform models are far less developed [23, 37, 38], 
because long, accurate simulations are particularly hard 
to achieve and because the parameter space is larger. The 
outcome and aftermath of NS-NS and NS-BH binary sim¬ 
ulations (as opposed to BH-BH ones) depend on several 
assumptions on currently underconstrained physics, such 
as the NS EOS, magnetic field configurations, and neu¬ 
trino emission. The relatively high mass ratios in NS- 
BH systems cause both analytical and numerical com¬ 
plications: the convergence of the PN approximation is 
slower than for NS-NS systems [39, 40], initial data are 
hard to construct [41] and the simulations must track 
very different dynamical time scales. In fact, BH-BH GW 
templates are commonly used in NS-BH merger searches. 
Tidal disruption affects both the GW and EM emission of 
NS-BH binaries. For all these reasons, better models can 


directly improve GW template banks, determine whether 
NS-BH mergers can power SGRBs and emit EM radia¬ 
tion in general (thus helping multimessenger searches), 
and improve event rate calculations for these systems. 
II. Simulations and model. At least two papers 
attempted a phenomenological description of the GWs 
emitted by NS-BH binaries. Lackey et al. [37] devel¬ 
oped an analytic representation of the IMR waveform 
calibrated to 134 numerical waveforms produced by the 
SACRA code [42] with the goal of assessing the measura¬ 
bility of the NS tidal deformability. The subset of sim¬ 
ulations for non-spinning BHs was used by Pannarale et 
al. [38] to build a frequency domain phenomenological 
waveform amplitude model which was, at heart, a “dis¬ 
tortion” of the PhenomG BH-BH model. This model 
relied on the fit of numerical-relativity results presented 
in [43] to compute the remnant torus mass Mb,torus- R 
paid particular attention to the accuracy at high frequen- 
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cies - where the EOS-related phenomenology takes place 
- and to the determination of the GW cutoff frequency. 
In units in which the total mass of the system is set to 
unity, let h{f) be the Fourier transform of the GW signal, 
and /Max the frequency at which f^h{f) is maximum. 
We define /cut (> /Max) as the frequency at which the 
dimensionless amplitude drops by one e-fold: 

c/cut^(/cut) — /Max^(/Max) ■ (1) 

Figure 1 displays /cut for a sample of BH spin parameters 
and for two piecewise polytropic FOS models (2H and 
B), chosen because they yield low- and high-compactness 
NSs, respectively (see [23]). This cutoff frequency is 
important to construct GW template banks for NS-BH 
searches: targeting NS-BH binaries with BH-BH tem¬ 
plates terminated at a frequency /term < /cut results 
in a signal-to-noise ratio loss; on the contrary, using 
/term > /cut may penalize the template by degrading its 
chi-square test performance, as it lacks matter effects. 

The companion paper [23] extends the work of [38] 
to NS-BH systems with a non-precessing, spinning BH, 
using the full set of 134 hybrid waveforms constructed 
in [37]. These are based on simulations in which the 
NS matter at zero temperature is modeled via piecewise 
polytropic EOSs that mimic nuclear-theory-based EOSs 
with a small number of parameters [44] . The binary mass 
ratio takes the values Q G {2, 3, 4, 5} and the BH di¬ 
mensionless spin parameter x G {—0.5, 0, 0.25, 0.5, 0.75}. 
This parameter space coverage allows our model to pro¬ 
duce the most accurate prediction of cutoff frequencies 
for NS-BH GW signals, with relative errors below 10%, 
well below the errors one would obtain using either BH- 
BH models or the NS-BH model of [37]. The new NS-BH 
GW amplitude model used here, and detailed in Ref. [23], 
is adapted to the three possible fates of the binary (see 
also Fig. 1): (1) nondisruptive: the GW frequency at the 
onset of tidal disruption /tide > /rd (where /rd is the 
BH remnant dominant ringdown frequency, calculated 
as in [45, 46], and /tide is the frequency at the onset of 
mass-shedding, determined as in [43]), and Mb,torus; com¬ 
puted as in [43], vanishes; (2) disruptive: /tide < /rd and 
Mb,torus > 0; or (3) mildly disruptive: either /tide < /rd 
and Mb, torus — 0, or /tide ^ /rd and Ibfb, torus ^ d- The 
BH ringdown does not contribute to the GW emission 
of disruptive mergers. Ringdown radiation appears in 
mildly disruptive mergers, and looks similar to BH-BH 
mergers in nondisruptive cases. Disruptive mergers thus 
differ the most from BH-BH mergers, precisely because 
tidal effects are maximal. In these events, /cut and the 
NS EOS have a strong link. To encompass all NS-BH bi¬ 
naries with a GW signal deviating from the BH-BH case, 
one must consider disruptive and mildly-disruptive merg¬ 
ers, i.e. binaries above the top dashed curves in Fig. 1 
must be discarded. The remaining set of binaries in¬ 
cludes the Mb,torus = 0 surface, and thus all possible EM 
sources (within the approximations of the model). 



FIG. 2. Boundaries separating disruptive mergers with 
/tide < /rd and Mb,torus > 0 (below each curve) from mildly 
disruptive and nondisruptive mergers (above each curve) for 
specific BH spin values indicated in the legend from top to 
bottom. Continuous (dashed gray) lines refer to EOS 2H (B). 

III. Predicting the fate of NS-BH mergers. The 

two dashed curves in the panels of Fig. 1 separate disrup¬ 
tive mergers (bottom-left region), non-disruptive mergers 
(top-right region), and mildly disruptive mergers (region 
in between the two lines). The green continuous line 
marks the locus of binaries for which Mb,torus goes to 
zero. Here we construct simple formulas to quickly de¬ 
termine the fate of a NS-BH coalescence. The contours 
in Fig. 1 that separate NS-BH binaries with a disrup¬ 
tive fate from those with a mildly disruptive or nondis¬ 
ruptive fate may be fitted in several ways as a func¬ 
tion of the binary physical parameters. We find it best 
to fit the critical mass ratio Qd below which mergers 
are classified as disruptive via a function of the form 
Qd = Qd(C)X)- Fitting in terms of the NS compact¬ 
ness C, rather than fixing an EOS and fitting in terms of 
the NS mass, allows us to use at the same time data pro¬ 
duced with the two extreme EOSs B and 2H. This is ev¬ 
ident in Fig. 2, where the boundaries between disruptive 
and mildly disruptive or nondisruptive binaries are shown 
in the CQ-plane for initial BH spin parameter values 
X G (-0.5,-0.25,0., 0.25,0.5, 0.75}: results for the 2H 
EOS (continuous) and the B EOS (dashed) overlay. This 
is due to the fact that the criterion to determine whether 
a binary is disruptive or not depends on Mb,torus; and 
that the state-of-the-art model for Mb,torus [43] depends 
only on C (i.e., it does not include higher-order, EOS- 
dependent effects). Quite independently of the EOS, the 
merger of an NS-BH binary will be disruptive whenever 
Q < Qd(C;X); where the threshold is well fitted by 
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Qd = ^ aijCx^ ; * + J < 3, (2) 

i,j—0 

with coefficients reported in the first row of Table I. The 
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TABLE I. Values of the coefficients of the fits in Eqs. (2)-(4). The number below each coefficient symbol must be multiplied 
by the corresponding power of ten in square brackets. The fijk’s are expressed in units of G = c = total mass = 1. 
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relative errors between the data in Fig. 2 and our fit are 
below ~ 4% [23]. Similarly, the relation 


/ 3 \ 
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Lv^=o / J 


with coefficient values listed in the second row of Table 
I, fits the boundary between nondisruptive and mildly 
disruptive mergers so that NS-BH systems with Q < 
Qnd(C, x) either disruptive or mildly disruptive. The 
maximum relative error between the data and the fit is 
^ 8%, and it is below 4.5% for 95% of the data points. 
IV. Cutoff frequency of disruptive mergers. Equa¬ 
tion (2) allows us to determine when a binary is disrup¬ 
tive according to our classification, i.e. the mass-shedding 
happens early enough during the evolution for the merger 
to produce a remnant disk mass and the GW emission 
deviates significantly from the BH-BH case at high fre¬ 
quencies. We now wish to provide a simple formula to 
compute the cutoff frequency of the GW amplitude /cut 
for disruptive mergers, as it carries information on the nu¬ 
clear EOS and it may be valuable in building better GW 
template banks. To this end we consider EOS 2H, gener¬ 
ate a set of 10^ random disruptive mergers, compute /cut 
for each NS-BH binary, and finally fit the resulting data. 
Disruptive mergers are selected as follows: we randomly 
sample parameters in the ranges M^^/Mq S [1.2,2.83], 
Q G [2,10], X G [—0.5, 0.75]; we verify whether the sam¬ 
pled point corresponds to a disruptive binary, i.e. whether 
/tide < /rd and Mb,torus > 0; we keep the point if it 
does; and we repeat the whole process until we have the 
desired 10^ points. While the maximum NS mass for the 
2H EOS is ^ 2.83M0, the maximum NS mass we obtain 
for the sample of disruptive NS-BH mergers is ~ 2.28Mq. 
The resulting mass interval Mns/Mq G [1.2,2.28] corre¬ 
sponds to a compactness interval of 0.117 < C < 0.221. 
We then fit the data set with the function 
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/cut= 51 i+j+k<8. (4) 

The resulting fijk values are reported in Table I. The 
relative errors between this fit and the data are typically 


below 1%: the relative error for 68%, 95%, and 99.7% of 
the points is 0.47%, 1.5%, and 4.9%, respectively. 

As a consistency check, we draw a separate sample of 
10^ disruptive mergers, and, for each binary, compute 
/cut and the relative error yielded by Eq. (4). This time 
we use EOS B, which has a maximum NS mass of ~ 2Mq. 
The compactness now ranges from 0.161 to 0.225. Quite 
remarkably, the maximum relative error is 2.2%; the rel¬ 
ative error of 97.6% of the points is below the percent 
level. This check implies that Eq. (4) is to a good approx¬ 
imation EOS-independent, at least within the parameter 
space where our model was calibrated and applied. 

V. Summary and discussion. In this Letter we used a 
recently developed semianalytical GW amplitude model 
for NS-BH mergers [23] to construct simple fits for: (i) 
the critical binary mass ratio Qu below which the merger 
is disruptive [Eq. (2)] and its GW emission deviates sig¬ 
nificantly from a BH-BH-like behavior; (ii) the critical 
binary mass ratio Qnd below which the merger is either 
disruptive or mildly disruptive [Eq. (3)]. This can be 
viewed as a necessary but not sufficient condition to gen¬ 
erate an EM counterpart, and may thus be used to deter¬ 
mine which binaries are plausible targets for multimes¬ 
senger searches targeting GWs and EM/neutrino emis¬ 
sion; (iii) the cutoff frequency /cut for disruptive mergers 
as a function of the initial binary parameters C, Q and 
X [Eq. (4)]. This can be used to maximize the recovered 
signal-to-noise ratio and chi-square test performance of 
BH-BH templates in NS-BH searches. The cutoff fre¬ 
quency can also be used to constrain the NS EOS [19- 
22, 47-49]. Our fit suggests that measurements of C, Q 
or X from the inspiral radiation (see e.g. [50, 51]) could 
improve the resulting constraints on the EOS. 

The non-negligible eccentricity of the initial data used 
for the NS-BH simulations underlying our model, the lim¬ 
ited duration and finite numerical resolution of the simu¬ 
lations, and the fitting errors all limit the accuracy of /cut 
to a few percent, and therefore introduce systematic er¬ 
rors. These errors are expected to increase as Q —t Qd, 
and when our fits are extrapolated beyond the region 
where the model and fits were tuned. More importantly, 
it is necessary to extend our model to processing bina¬ 
ries [52] . We plan to address these issues in future work. 
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